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Abstract 

We study the criticality of a Potts interface by introducing a froth model 
which, unUke its SOS Ising counterpart, incorporates bubbles of different 
phases. The interface is fractal at the phase transition of a pure system. 
However, a position space approximation suggests that the probability of 
loop formation vanishes marginally at a transition dominated by strong ran- 
dom bond disorder. This implies a linear critical interface, and provides a 
mechanism for the conjectured equivalence of critical random Potts and Ising 
models. 
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The effect of quenched impurities on phase transitions is important and quite fascinat- 
ing. A simple "Harris criterion" [|I| indicates that critical behavior is modified by (bond) 
randomness in systems with a positive heat capacity exponent. Impurities can also change 
first order transitions to second order 0, to the extent that in two dimensions there are 
no discontinuous phase transitions [Q. There is a growing body of numerical |Q and exper- 
imental P] evidence that, at least in some situations, the asymptotic criticality is similar 
to the random bond Ising model, irrespective of the underlying symmtery. Here we provide 
some justification for this observation based on an unexpected universality of the critical 
interface in the presence of strong bond randomness. 

Interfaces are key to second order phase transitions; the interfacial free energy vanishes 
at criticality with the Widom exponent of fi. The behavior of Ising interfaces has been ex- 
tensively studied by a solid-on-solid (SOS) model which simplifies numerical and theoretical 
analysis [0. Potts models also have a potential continuous transition with broken discrete 
symmetry 0. Potts interfaces are harder to study as they are complicated branching objects 
which become fractal at criticality. We introduce several approximations to such interfaces 
in 1 + 1 dimensions, eventually arriving at a model that is amenable to a position space renor- 
malization group (RG) treatment; exact on a hierarchical lattice. The simplified model also 
allows us to examine interfaces of random bond Potts models. This leads to an intriguing 
generalization of directed polymers in random media p|, with potential applications beyond 
those discussed here. 

Consider the interface between two distinct phases of an ordered q state Potts model. At 
very low temperatures the interface is a weakly fluctuating surface, much as in the Ising case. 
However, on approaching the critical point, bubbles of any of the other q — 2 phases may 
appear at the interface, in addition to islands in the bulk phases. As in SOS models of Ising 
interfaces, we shall ignore isolated islands and overhangs. The resulting froth is a collection 
of bubbles, each bounded by two SOS surfaces. Even this simplified model is too hard to 
analyze: The allowed configurations in 1+1 dimensions are a subset of those encountered 
in directed percolation (no dangling branches), with weights depending on q. To make the 
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problem tractable, we confine the interface to the bonds of a diamond hierarchical lattice 
(DHL) (inset of Fig.(l)). In one iteration on this lattice the interface can cross either of two 
branches, or create a bubble by going through both. The bubble has a "fugacity" of q — 2, 
the number of possible intermediate phases. This procedure is repeated iteratively creating 
the possibility of loops within loops ad infinitum. Not all configurations of the original 
froth model (those with multiplicities that are not powers of g — 2) are included within 
this scheme. It is possible to construct more complicated models and RG schemes without 
these deficiencies; for example by considering the diamond plus diagonal hierarchical lattice 
(DDHL) in Fig.(l). We have checked that the qualitative results are unchanged by the 
choice of RG scheme, and will thus focus on the simpler DHL where the recursion relation 
for the interface partition function is 

Z„+i = 2Z^ + (g-2)Z^ . (1) 

Here Zi = exp(— 2/5J) is the Boltzmann weight of one broken bond and (3 = l/{kBT). 

The above recursion relation has a stable fixed point aX Z* = 0. For initial values 
flowing to this point Zn oc exp[— L/s], where L = 2" is the length of the lattice after n 
iterations and fg is the interfacial free energy. There is another stable sink with Z* ^ oo 
where Z„ oc exp[— L^/^]. This phase is a dense foam of bubbles where the interface analogy 
breaks down. These two phases are separated by an unstable fixed point at a finite Z*{q). 
Clearly any finite fixed point for Z corresponds to /s = = 0! A similar fixed point 
mechanism for the vanishing of the surface energy is present in wetting phenomena 0. 
Denoting t = Zi - Z*{q), we find t' = 2y^iH with y{q) = 2 + ln(l - Z*(g))/ln2. Since the 
interface free energy satisfies the homogeneity condition fs{t) = b^^fs{by^'^h), it vanishes as 
with yu(g) = l/y{q). The exponents obtained from the DHL and DDHL RG schemes 
are compared with the exactly known values for the Potts model in d = 2 in table I. The 
exponent is (accidentally) exact for g = 2 on the DHL, and shows the correct trend on 
increasing q. However, this approach does not show the expected change to first order 
transitions at g = 4 0. The critical interface is fractal, its mass (number of occupied bonds) 
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growing as M = d\nZ/d\nZi oc L'^f with df = y{q). The probabihty of forming a loop is 
independent of scale and given by P = — 1. 

As shown by Hui and Berker [0, randomness in the couplings of a Potts model may 
change first order transitions to second order. In particular, there can be no discontinuous 
symmetry breaking in two dimensions , and all disordered Potts models have continuous 
phase transitions. This opens up the possibility of an infinite set of new universality classes. 
A numerical study 0] of the 8-state Potts model has confirmed the second order nature of 
the transition. Intriguingly, based on their simulations, the authors of ref. conjecture 
that all these new universality classes are in fact similar to the random bond Ising (g = 2) 
model. There is also some experimental support for this conjecture from the effects of 
oxygen impurities on the phase transition of (2x2) — 2H on Ni(lll). The pure system is 
a realization of the g = 4 Potts model, while in the presence of impurities the asymptotic 
behavior appears to be Ising-like . This contradicts other evidence based on RG schemes 
p!0| [pJ]] that, at least for weak disorder, the critical exponents do depend on q. In fact 



the latter RG scheme |rT| was partly motivated by an earlier work suggesting that the 
(nonuniversal) critical line of the pure Baxter model shows Ising behavior in the presence of 
randomness. 

To gain further insight into random bond criticality, we examine the singular behavior 
of the interfacial energy in the framework developed earlier. The initial weights {Zi{i)} 
for the bonds of the DHL are chosen randomly from a probability distribution Vi{Z). The 
recursion relation for a specific realization of bonds, 

Z„+i= Z„(1)Z„(2) + Z„(3)Z„(4) 

+ (g-2)Z„(l)Z„(2)Z„(3)Z„(4), (2) 

can be used to construct a functional recursion relation for Vn{Z). We use a binary initial 
distribution of energies with a fraction p of positive bonds Ji, and 1 —p of positive or negative 
bonds J2, i.e. 

Vi{Z) =p6(Z- e-^f^A + (1 - p)6 (Z - e-^^^A . (3) 



Under iteration, the simple fixed points of tlie pure system are replaced by stable distri- 
butions: The analog of Z* = describes directed paths on the hierarchical lattice and is 
governed by the distribution discussed in ref. [|13|. This distribution for InZ is characterized 
by a mean that shifts as —fsL <^ and a width that grows as L'^ with uj ^ 0.30 (compared 
to the exact value of = 1/3 on a two dimensional lattice M). In the dense foam phase with 



\nZ = —fbL"^ ^ 0, the central limit theorem should apply and we expect fluctuations of 
In Z to grow as L. An argument similar to the Harris criterion [|l| shows that the relevance 
of randomness at the critical fixed point is determined by the sign of y{q) — 1. Thus any 
randomness is expected to modify the interfacial criticality of Potts models with q > 2. 
By analogy with the other two limiting distributions, we expect the critical point to flow 



towards a third zero temperature stable distribution with InZ = 0, and S\nZ oc L^". 

These expectations are consistent with numerical iterations of the recursion relation 
starting from a large initial ensemble of {Zi{i)} The RG flows are towards zero tem- 
perature, leading to the phase diagrams indicated in Fig.(l) for J2 = —Ji (and J2 = +Ji/2) 
for q = 3 and 4. The choice of an initial distribution with a fraction of negative bonds 
ensures that the phase boundary extends to zero temperature Since the eventual fixed 



distributions are at zero temperature, criticality can be examined by directly looking at the 
recursion relations at T = 0. With this enormous simplification, we have iterated energies, 

En+i= min{E„(l) + En{2), E„(3) + E„(4), 

E„(l) + E„(2) + E„(3) + E„(4)}. (4) 

The last term is of course absent for q = 2. We could iterate exactly the probability 
distribution for bond energies up to n = 7. The Monte Carlo iteration was typically extended 
up to n = 20 for our determinations of the T = critical properties. It is not immediately 
apparent from Fig.(l) that the phase boundary for an initial mixture of ferromagnetic bonds 
is also be governed by a T = fixed point that has a mixture of positive and negative bonds. 
Although in a regular RG scheme positive bonds only generate positive bonds, it can be 
checked easily that for all temperatures intermediate between the extreme critical points. 



the recursion of eq.(2) leads to a mixture of positive and negative bonds. Thus the critical 
behavior is likely to be the same for both cases. 

A striking feature of the T = recursion relation is that it is independent of q for q 2. 
Thus the critical behavior for the vanishing of fs should not depend on q. This expectation 
goes beyond the approximations of the model and the hierarchical lattice. Consider the 
configurations contributing to a low temperature expansion of the interface free energy 
of the Potts model on any lattice. The index q appears only in entropic factors giving 
the multiplicity of possible colorings. If, as is usually the case in random systems, the 
scaling properties of the interface are controlled by a zero temperature distribution, these 
properties will be independent of q\ Of course, the SOS configurations allowed for g = 2 
are very different from the froth that appears for q > 2, and thus in principle, we expect 
two different universality classes. Strictly speaking, since four colors are needed to cover 
an arbitrary "map" in two dimensions, there are restrictions on configurations allowed for 
g = 3. If relevant, this leads to a third potential random interface behavior in d = 2. Indeed 
the recursion relations of the DDHL at T = discriminate between g = 2, g = 3, and g > 4. 

We focused our studies of critical behavior at T = on the mixture of bonds of strengths 
+ J and — J. The critical concentration in the presence of bubbles was identified by trial and 
error as pc = 0.883 ± 0.001. At the critical point, the mean and variance of energy should 



scale as E{L) = L^" and var[i?(L)] = L'^^", with O^. = ujc if there is only one energy scale. The 
fits at Pc are consistent with this expectation, giving O^. = 0.34 ± 0.05 and Uc = 0.31 ± 0.08. 
Rather surprisingly, Uc is very close to the value of u for directed paths with no loops. 
Fig. (2) shows that the data for the energy on approaching pc from above can be collapsed 
by a finite-size scaling form, 

Emp>pc) = A\pc-prL + BL'^% (5) 

with fi = 0.91 ± 0.04. In the absence of loops for g = 2, /i = 1 exactly, since the mean value 
of the final energy is simply linear in the mean value for individual bonds. There is also a 



subleading correction to the energy from fluctuations that scales as L^" ||T6|. 



To better understand the closeness of exponents in the presence and absence of loops, 
we looked directly at the fractal structure of the critical interface. This is achieved by 
examining the probability P„ that a loop forms at the n}^ iteration. As indicated in Fig. (3), 
Pn decays exponentially to zero for p > pc, while the decay slows down on approaching pc- 
The data at pc can be fitted to a decay as c/n, with c = 0.060 ± 0.003. This fit suggests the 
phenomenological differential recursion relations 

The first equation reproduces P{n) = c/n at criticality; the second is just the behavior of 
the mean (free) energy (or mass of the cluster). The fixed point at /<j = P = describes 
a critical cluster that is asymptotically linear, with loops appearing predominantly at short 
length scales. Linearizing the second equation gives y = I. However, the marginality of P 
at criticality leads to logarithmic corrections to various scaling quantities. For example, the 
mass of the critical cluster grows as L(lnL)'^. Similarly, the interfacial free energy vanishes as 
|p— Pc|| lii(|p— PcDI'^- Fig- (2) shows that a data collapse is also possible using such logarithmic 
corrections. The best fit is achieved for c = 0.5±0.6, compared to c = 0.060±0.003 obtained 
directly from Fig. (3). 

The asymptotic linearity of the critical cluster extends beyond the simple example of 
the DHL, and was also checked for the DDHL. In the latter, an interface configuration that 
covers both branches, as well as the central diagonal, is allowed only for g > 3. Thus this 
lattice supports three types of recursion relations at T = 0, corresponding to g = 2, 3, 
and > 4. In spite of these differences, within our numerical accuracy, we could not detect 
any significant changes in the exponents from the simpler DHL. As indicated in table I, 
even in the absence of loops {q = 2), the critical cluster in the pure system is a fractal. 
This is because it takes advantage of the diagonal bond and is no longer simply directed. 
The critical cluster of the random system does not take advantage of the diagonal bond, 
or the possibility of loops, staying asymptotically linear. This suggests that the marginal 
irrelevance of operators that may complexify the structure (such as P in eq.(P)) is in fact 



quite generic. 

Although we have presented the resuhs in the context of Potts models, they are probably 
more generally applicable to systems of discrete symmetry. For example, it has been sug- 
gested that regions of the Baxter model with diverging heat capacity also exhibit Ising like 
behavior in the presence of random bonds [|T^ . It is also tempting to generalize the con- 



clusions to higher dimensions: Strong enough disorder may result in a continuous transition 
with the interfacial criticality governed by a zero temperature fixed point. If the appar- 
ent irrelevance of bubbles can be generalized from hierarchical lattices, the exponent n will 
be the same as in the Ising model. How can we reconcile this apparent super-universality 
of random bond criticality, with earlier results |0 which do indicate exponents that 
depend on g? The latter calculations were performed for weak disorder, and lead to a finite 
temperature fixed distribution. Is it possible that stronger disorder leads to different behav- 
ior, dominated by zero temperature fixed points? This scenario is precisely what is observed 
in a recent position space RG of a three state random bond model |]18[. Of course, another 
possibility is that our approximate interface model does not fully capture the physics of the 
system. Clearly, further investigations are desirable. 

In summary, we have introduced a simple position space approximation for studying the 
interfacial properties of Potts models in d = 2. In the pure case, the critical interface is a 
fractal froth, and interfacial tension vanishes with an exponent fj^{q). In a random medium, 
within our numerical accuracy, we find that the critical interface is asymptotically linear, 
and the interfacial free energy vanishes linearly with logarithmic corrections, i.e. /i = 1 
independent of q. We argue that, if governed by a zero temperature fixed point, fi{q) should 
be independent of q for all sufficiently large q. The above model also provides the simplest 
generalization of directed paths in random media p] to ramified objects, with potential 
applications to fracture cracks, lightning patterns, etc. 
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TABLES 

TABLE I. The exponent jLt(g) for the vanishing of the interfacial free energy, obtained from the 
DHL and the DDHL, compared to the exact values in d = 2. 



q DHL DDHL Exact 

1 L634 1.429 4/3 

2 1 0.847 1 

3 0.886 0.749 5/6 

4 0.830 0.700 2/3 
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FIGURES 



FIG. 1. Phase boundaries for g = 3 and 4, with mixtures of positive and negative bonds (solid, 
dashed), and positive bonds (dotted, dashed dotted). Arrows indicate RG flow towards the T = 
fixed point. The solid lines in the inset figure indicate the cell replacing each bond at every stage 
of the DHL construction. The central (dashed) diagonal is also present in the DDHL. 



FIG. 2. Data collapse of E{L,p > Pcj/L'^", against {p - pcY L'^'^-'^'^'i (left group), or 
\p — Pell ln(|p — Pc\)\'^L^^'~''^''^ (right group) in the abscissa. The scales are logarithmic. 



FIG. 3. The iterated probability for loop formation P„. The decay is fitted to c/n (asterisks 
and dashed line) at Pc- Data for p > Pc (squares) can be fitted to an exponential. 
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